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Some brane models rely on a generalization of the Melvin magnetic universe including a complex 
scalar field among the sources. We argue that the geometric interpretation of Kip.S.Thorne of 
this geometry restricts the kind of potential a complex scalar field can display to keep the same 
asymptotic behavior. While a finite energy is not obtained for a Mexican hat potential in this 
interpretation, this is the case for a potential displaying a broken phase and an unbroken one. 
We use for technical simplicity and illustrative purposes an ad hoc potential which however shares 
some features with those obtained in some supergravity models. We construct a sixth dimensional 
cylindrically symmetric solution which has two asymptotic regions: the Melvin-like metric on one 
' side and a flat space displaying a conical singularity on the other. The causal structure of the 

, configuration is discussed. Unfortunately, gravity is not localized on the brane. 
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I. INTRODUCTION 

> 

Among the most important characteristics of cosmic strings is the existence of a symmetry axis and the concentration 
of energy around this axis Ij. Taking gravity into account, the existence of a symmetry axis implies cylindrical 
symmetry for the metric as well. The static cylindrically symmetric solutions of Einstein equations in vacuum in 4Z? 

• are of the Kasner type: they are parameterized by three constants obeying two constraints. The vanishing of the 
' energy-momentum tensor in the asymptotic region implies that the geometry must approach a Kasner line there. As 
: the energy momentum tensor corresponding to these axial configurations implies the invariance of the metric under 
' boosts^ne is left with only two Kasner geometries : a flat space presenting a conical singularity and a Melvin like 

> , space 01 ■ The first case leads to the well known cosmic strings. The second solution, written in a particular system 

• of coordinates, displays circles of decreasing circumferences for increasing "radii" p. This feature has been analyzed 
by Kip.S.Thorne The interpretation is that p = 00 is the point at infinity on the symmetry axis. 

The Melvin solution has high dimensional generalizations which can be used to build brane models 0, 0| . In recent 
^ works, complex scalar fields have been incorporated into the picture 1^. In this article we address the same question. 
However, we use the Kip S.Thorne interpretation to fix the boundary conditions; the coordinate p being the point at 
infinity on the symmetry axis, the angular coordinate is not well defined there, just as for the polar coordinates at 
the origin on the plane. A cylindrically symmetric complex field must thus vanish at that point. We obtain that no 
static cylindrical solution can be obtained with the usual Higgs potential. We then exhibit a toy model for which this 
can be achieved and discuss its characteristics. 

This paper is organized as follows. In the second section we review the geometric interpretation of The Melvin 
magnetic solution in four dimensions. The third section incorporates the scalar field, set the boundary conditions 
and displays the numerical solution obtained. First, an Abelian-Higgs Lagrangian is coupled to the Einstein-Hilbert 
one. Looking for an axially symmetric configuration which displays the second special Kasner line element far from 
the source , the preceding section imposes the vanishing of the vector and the scalar field as the coordinate p goes 
to infinity. This results in a divergent inert mass per unit length if the v.e.v of the Higgs field does not vanish. 
If it does, one simply recovers the Melvin solution. We then construct, for illustrative purposes, a potential for 
which the behavior of the scalar is non trivial. This potential has two minima, one of them being at zero. This has 
similarities with some potentials obtained in some supergravity models and some non commutative models |^. 
The constructed solution interpolates between them. The minimum corresponding to a vanishing value of the field is 
attained in the Kasner-like asymptotic region while the non vanishing value corresponds to a fiat space presenting an 
angular deficit. The classical trajectories of neutral particles in this geometry are analyzed in section four. We show 
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that the trajectories of massive particles in this geometry are bounded. We also study massless particle trajectories 
and discuss the causal structure of the solution. An appendix is devoted to the way the numerical approximation has 
been computed. 

The coupling of scalar fields to gravity leads to many classical solutions [3. ^i IT^ ll^ ITsl H^ . The 
introduction of a scalar field among the sources leading to a Melvin type universe has been considered before [T^ . 
There are three main differences between our work and the previous papers. Firstly, the scalar field considered here 
is not a dilaton so that its coupling to gravity is ordinary. The involved potentials are not the same. Secondly, the 
scalar field here is complex, contrary to JJj where it is real. The vanishing of the scalar field on the symmetry axis 
is not required when it is real. On the contrary, this becomes mandatory when it is complex, just like for the Higgs 
field on the cosmic string core. But, the more important difference is that here we have a configuration with two 
asymptotic regions. Like in f^, we do not have a delta function for the brane. 



II. THE MELVIN SOLUTION IN AD. 



In this section we review the geometric interpretation of the Melvin solution. To keep the discussion as simple as 
possible, we actually analyze its asymptotic limit in four dimensions. The conclusions are however the same in the 
presence of extra dimensions and the addition of matter. 

We will study a system consisting of a self gravitating Maxwell system in d dimensions. The classical field equations 
are derived from the action 



S= d'' 



1 p pab , R 



4 167rG 



(1) 



The solution which generalizes the Melvin universe in d dimensions is given by the following expressions of the metric 
and the Maxwell tensor: 



ds^ ^ + (^rj^^dx^dx''-dp^)-p^l + ^^ 



2\ -2 

P 



F - Bo 1^1 + pdpAd(b , (2) 

a and Bq being dependent constants. This solution has been used in the brane world scenarios. One of its important 
properties is that the brane can have positive tension and the closure of the bulk provides a singularity-free boundary 
condition for solutions that contain black holes and gravitational waves on the brane 0, . 

The characteristic of this metric on which we will put the emphasis is the fact that the circumference of circles 
obtained by letting only </) non fixed tend to zero as the coordinate p goes to infinity. The Lagrangian from which the 
solution given in Eq.Q follows did not contain any scalar field. The question we address in this paper is: if one adds 
a scalar field to the picture, which kind of potential allows the same behavior for the metric? We will argue that a 
geometric interpretation of this property of the metric gives an important restriction. 

The features observed bv Kip.S.Thorne for the Melvin solution IS., .19j are already present in a vacuum solution: 
the Kasner line element which is obtained by taking a = in Eq.Q. The simplest example where this behavior 
can be analyzed is the sphere. Introducing the polar stereographic coordinates (r, 9), the metric of a sphere of radius 
a reads 

16cr^ lgfj4j,2 
(4(t2 + r2)2 (4cr2 ^ j,2)2 ^ ' 

For large values of the coordinate r, the coefficient ggg becomes a decreasing function . If one interprets r as the 
radius, a circle of infinite radius turns out to be of null length. This is obvious since r = oo corresponds to the point 
at infinity on the plane which is mapped into the north pole by the stereographic projection; r = oo is just the north 
pole. When the coordinate r vanishes, one has another circle displaying a vanishing circumference : the south pole. 
The two are on the symmetry axis. Introducing the variable 

r* = 2CTarctan(r/2(T) , the metric reads ds^ — dr^ + s\t? dO^ . (4) 



The relation between r, and r is bijective provided that € [0,7rcr]. The points located on the symmetry axis once 
again are those for which the coefficient ggg vanishes. 
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III. THE 6D EXTENSION WITH A COMPLEX SCALAR FIELD. 



Before proceeding, let us point out that the tangential Maxwell vector field in the Melvin solution vanishes when 
p goes to infinity, in accord with the geometric interpretation. We now wish to include a scalar field in the picture, 
in the presence of extra dimensions. We choose a sixth dimensional model essentially because in this case one can 
naturally obtain chiral fermions. 

Let us first consider a scalar displaying a Higgs potential; the matter action is then 



S = 



-9 



2 ^ 4^ ^ 



(5) 



The U{1) charge e is embodied in the covariant derivative 1?^$ = 3^$ — leA^^. For a static cylindrically symmetric 
configuration, the ansatz can be given the form 



ds^ = (3^{p)r]^^dx^'dx'' - -/^{p)dp^ - a^{p)d(jp- where /i = 0, • • • 3 
<t> = v!{p)e"^ and A^^-{l-p{p)) . 



(6) 



The cosmic string solution has been extensively studied in the literature. In that configuration, the smoothness of the 
eometry on the symmetry axis is guaranteed by going to the gauge 7(p) = 1 and imposing the boundary conditions 



(7) 

(8) 



a(0) = 0,a'(0) = 1 , 
while the matter fields are non singular on the core provided that 

/(O) = , p(0) = 1 . 

The vanishing of the energy density in the asymptotic region implies 

/(oo) = l , p(oo)=0 



(9) 



What happens if we want the metric to display the same asymptotic behavior than in Eq.(|2Jl when the coordinate 
p goes to infinity? In the previous section, we argued that p = cx3 is the point at infinity on the symmetry axis. To 
have a regular cylindrically symmetric configuration, the Higgs and the tangential vector field must vanish there: 



/(oo) = and p(oo) = 1 
Extracting the expression of the integrand of the inertial mass from Eq. jSJ one has 



\9\ 



2\2 



(10) 



(11) 



In the asymptotic region (i.e. p co) one has ^/— 5 p'^^\ the volume element is not bounded. The first three 
terms decrease in the asymptotic region provided that /(p) and p(p) approach constants there; this is already satisfied 
by Ea. pOfl . The contribution of the Higgs potential in this part of the space is reduced to the integral of p'^^^Xv'^. 
This has a chance to converge only when v = Q. Then, the parameterization given in Eq.@ does not apply; one can 
however define a dimensionless function associated to the Higgs field by using the Newton constant. Doing this, we 
obtained a vanishing scalar field for any value of the parameters. 

Physically this can be understood as follows. Forcing the scalar field to go from zero to zero as p goes from zero to 
infinity, one obtains that it vanishes identically since there is no source. Such a source would be for example a local 
maximum of the potential but as the vacuum expectation value vanishes, such a maximum does not exist. In fact, 
one recovers the Melvin universe. 

Is it possible to construct a solution with a non vanishing scalar field? To do this we need a potential which vanishes 
with the scalar field so that the minimum is attained at spatial infinity. We also need a local maximum which will 
correspond to a source. These conditions are for example satisfied by the gauge invariant potential 



(12) 



The maximum is attained at = ±u\/v2— 1 while there are three minima , at ^ = 0,±'!;. The U(l) symmetry is 
broken spontaneously in the last two vacua and preserved in the first one. We disregard the renormalizability since 
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we are interested only in classical solutions. This potential, although purely ad hoc, shares with the one appearing in 

Q 

V{(j)) = 2e^*(#)f-i [2(p + - 500] (13) 

the fact that it is the product of an exponential and a polynomial. The difference is the fact that in Ea. (|i;-i|l . one needs 
to have p > 1 to have the vanishing value of the field as a vacuum but then there is no other vacuum. In , it was 
argued that the potential of Eq. H13|l could be seen as inspired from some supergravity model, with a particular choice 
of the Kahler structure. Let us also mention that in models in which one works with non commutative spaces such as 
A4 X Mn where Ai denotes the Minkowski space and M„ the set oinxn matrices, one also obtains potentials displaying 
symmetric vacua. The attitude adopted here is like the one concerning analytical solutions for self gravitating domain 
walls p^. I23l l2^. One knows that a potential which is a cosine of a scalar field achieves the desired goal, although 
it is not renormalizable. In the same way, one introduces ad hoc potentials for quintessence (25- 26, 27] . Some of 
them are negative powers of a scalar field and so lead to non renormalizable theories. Our only aim is to show that 
solutions with finite energy exist for specific potentials. Moreover, some potentials displaying symmetric vacua have 
been used as candidates for dark matter 28, 29] . 

We now wish to construct a solution which interpolates between two vacua, say |$| = w and |<i>| = 0. Our previous 
discussion tells us that the region where the field goes to its unbroken phase can not be Melvin-like. As we simply 
require cylindrical symmetry, we can choose that asymptotic region to be like the far region of a cosmic string. 

The Einstein equations will be written in the form 

Rl = -8^G - ^t) (14) 



where the factor 1/4 comes from the dimension of the space time. In sixth dimensions, asking our action to be a pure 
number means the self coupling A and the gauge coupling e are dimensionfuU. We thus can write our equations in 
terms of the dimensionless parameters 

^ = Gt;2 , 1^ = X^G^^ , T^e^v and a = w'^v^ . (15) 

Note that here the dimension of v is an inverse length square. We will use the dimensionless length x given by 

p = Lx where L = ^^=-ij-^-L (16) 

and the dimensionless functions 

a{p) = ^.[l^A{x) , m=B{x) , ^{p)=i{x) , f{p) = F{x) , p{p) = P{x) . (17) 

V T 

The field equations read 

.-<..5,.).^Wf ,.,(1 - .^W,^ . ^iigM + 3^ - i» - . Bi.) = 

P''^(x) 

+6^^ + A"{x)=0 , (20) 



T 



FHxm^)Pix) ^ + _ + p"(,) . o , (21) 



27r^5/27I7 ^ " ^ ^ ^ A{x) B{x) 7(2;) 
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FIG. 1: The scalar and vector fields are plotted in terms of the coordinate x. 



5i;e-=<"F,.)f (.)(1 + (-4 + ,)FHr) + (3 - 2a)F'i.) + - j-^Silg^pl 

A{x) B{x) 7(2;) ^ ^ ^ ^ 

We now relax the assumption that the coordinate p goes from zero to infinity but rather take it to go from minus 
to plus infinity. This is like parameterizing every point of the sphere not by 6, <p but by (j) s^nd its distance from the 
equator on a meridian. The upper hemisphere would have positive p while the lower would correspond to negative 
values of that length coordinate. In our case, the region where p — oo will correspond to a cosmic string-like 
geometry while p — > oo will be associated with the Melvin-like behavior. 

So, the boundary conditions are that when p — oo, one is in the far region of the cosmic string solution: 

A{x) - x , B{x) - 1 , 7-1 (23) 
while for p oo one enters the asymptotic region of the Melvin solution: 

Aix) - - , B{x) - 7(^) - f 4) ■ (24) 

X \0' J 

To give an illustration, we need to solve the above set of coupled non linear differential equations. This has been 
done for the simplest choice of the parameters: p = v = T = ij = \] the details concerning the numerical treatment 
are given in the appendix. The behavior of the different fields has been given in the pictures FIG.l, FiG.2, and FiG.3. 
Roughly speaking, the region p = is where the transition between the two regions takes place; it is also the place 
where the energy is concentrated. 

The trapping of gravity is ensured when the condition 

„4 J 5 00 



dx*dx^5"Vl5l < oo (25) 

is satisfied. From the asymptotic behavior of the metric, one sees this is not true for the solution we constructed. Let 
also remark that a change of coordinate can now be made so that for example 7 = 1. 



IV. THE CLASSICAL TRAJECTORIES. 



Let us first consider the case of massive particles. The metric specified in Eq.© has cyclic coordinates ; its geodesies 
are by way of consequence characterized by constants of motion. Introducing the proper time per unit mass r, the 
energy-momentum relation reads 

dpV f kl k2 \ I 

d^J -['^'^^Wip))^) ■ ^^^^ 

A physical motion is characterized by a real velocity. The asymptotic behavior of the metric displayed in Eq. H23I24|) 
shows that the trajectory of a massive particle never attains the point at infinity on the Melvin branch; however it 
has access to the region at infinity on the string branch. 
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FIG. 2: Prom top to bottom, the components A{x),B{x),^{x) of the metric 
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FIG. 3: The energy density t{x) 



The causal structure is found by analyzing particular null geodesies. The bounded null coordinates in the new 
background are given by u = c** or v = c** with 



arctan[t =f (^{p)] 



where 



a{p)= r 
Jo 



(27) 



(28) 



In FIG. 4, we have drawn the Penrose-Carter diagram of the solution. 

The embedding of the two dimensional metric containing only t and p in the Euclidean space can be realized by 
the surface of revolution 



Z{r) = dyVcHy) - 1 , c{y) = j{a-\y))-{a-\y)) . 

The limiting behavior of the metric shows that in the region p — > — oo, 

Z{r) ~ c"*r ; 

no restriction is imposed on r and the surface is a cone. On the contrary, as /> — > oo. 



rr r / \ 16/3 ] 

Zir^.jjy (-) -I 



(29) 



(30) 



(31) 



there is a maximal circumference. 
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FIG. 4: The causal structure of the solution is given. The time-like infinities are /'''(it = v = 7r/2) and I~ {u — v — —n/2). As 
r can change sign, there are two space infinities : 7^(— 7r/2, 7r/2) and 7^ (7r/2, — 7r/2). The curves which begin at I~ and end 
at correspond to fixed values of r while the others correspond to fixed t. 

V. CONCLUSIONS. 

We have constructed a sixth dimensional cylindrically symmetric self gravitating configuration. It has two asymp- 
totic regions: one corresponding to a flat space time with a deficit angle and the other to the second special Kasner 
line element. The part of the geometry which displays a deficit angle can be realized as a cone with a smoothed apex. 
The second special Kasner geometry, on the other side, can be seen as a tube with a decreasing radius. Gluing the 
two, one obtains something close to a funnel. The causal structure of this geometry was studied. It should be stressed 
that such a configuration, with two topologically different boundary regions, is not possible with static spherically 
symmetric configurations; the Birkhoff theorem forbids this. 

The potential considered is non rcnormalizable but our discussion shows that it is one of the simplest which allows 
boundary conditions compatible with the geometric interpretation of the second special Kasner geometry. 

Like in |32|. we have built a configuration which has two different asymptotic regions. In our case we have a space 
with a conical singularity on one side and a Mclvin-like solution on the other, while in the preceding one there are two 
Ads5 space times glued together along a three brane. The fact that our solution is not asymptotically flat and the non 
localization of the four dimensional graviton is similar to j33j . Among the priorities which should be addresses if a 
more realistic model is built along these lines is the construction of realistic Abelian and non Abelian four dimensional 
models [3II33. 

What we have learned in this work is basically that if one wants to include a complex scalar field possessing a 
winding number on a Melvin-like solution, one needs a particular kind of potential. Our solution docs not trap 
gravity. Nevertheless, the model may still have some phenomenological interest. Although we have not made here the 
appropriate analytic computations, one can not rule out at this stage the possibility of having a quasi localized four 
dimensional graviton on the brane like in |36|. In this model, it was shown that Newton's law of gravity was valid 
only between two length scales fixed by the theory. 
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APPENDIX A: NUMERICAL CONSIDERATIONS. 

Our numerical approximation relies on a symbolic approximation of the fields which can then be improved by a 
relaxation method. The function F{x) goes from 1 to when the argument goes from — oo to oo. The ansatz is taken 
to be 

J^(^)-^(l-tanh(/gx)) f + ff . (Al) 

In the same way, one has 

P{x)^hl + tanh{plx))^f±4^ . (A2) 
2 1 + P2X'' 

The asymptotic behavior of the metric displayed in Eas. H23l24|) is taken onto account by the following functions: 
A{x) = ^{1 - tanh {alx)){al + alx'^)^/'^ + ^{1 + tanh {alx)){al + alx'^y^/'^ , 
B{x) = 1(1- tanh (6^.))|^t^ + 1(1 + tanh (6l.))(6l + 6^a;^)V3 , 

7(x) = 1(1- tanh (,g„^a:))|^i^ + 1(1 + tanh (,g|:.))(5^+ 55^)1/3 . (A3) 

These parameters are not all independent, due to the fact one constant(a) drives the asymptotic behavior of the 
functions a{p), f3(p),j{p) and P{p) simultaneously (see Eq.|(21)). This has been taken into account. For simplicity, we 
introduce 0:5 = -^05. 

For the true solution, all the right members of the equations given in Eg. ljlSll . • Ea. (|22|l , which we denote 
ODEi{x), ■ ■ ■ jODE^^x), vanish. To obtain an initial approximation for the relaxation method, the idea is to look 
for the values of the coefficients /o, • • ■ ,55 for which the integral 

/oo 
dx (ODEfix) + ODEl{x) + ODElix) + ODEI{x) + ODEI{x)) (A4) 
-00 

is minimal. 
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Rather than computing the integral, we approximated the surface to which it corresponds by a sum of rectangles. 
The values of the parameters we found are given below. Plotting the functions ODE].{x), one finds an error of the 
order 10~^ on the entire real axis. 

fo = 0.29525844042277477 , /i = 0.9658023897961495 , /2 = 0.21489954200347053 , po = 0.24041269454525904 , 
pi = 0.9931697892207005 , p2 = 0.4970289653744458 , oq = 0.8317552649976594 , ai = 1.3064942111721376 , 

a2 = -0.011436440955002887 , as = 0.8012738041970792 , 0,4 = 0.6275485649917462 , Q5 = 0.019766866475879406 , 
60 = 0.823008493904919 , h = 0.9933291170842984 , 62 = 0.42807374544606963 , 63 = 0.7966722551728198 , 
64 = 1.2646707186251818 , go = 0.7842562408108704 , gi = -0.2977638255573934 , 52 = -0.0349623024909614 , 
c/3 = 0.8326686788126545,54 = 0.14292786831126408 . (A5) 



